The Hg molecular ion in a magnetic field: linear parallel 

configuration 

A. V. Turbiner,* N. L. Guevara,^ and J. C. Lopez Vieyra-f 
Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 70-543, 04510 Mexico, D.F., Mexico 
(Dated: February 2, 2008) 

Abstract 

A first detailed study of the ground state of the HJ^ molecular ion in linear configuration, parallel 
to a magnetic field direction, and its low-lying S,n,A states is carried out for magnetic fields 
B = ^ — 4.414 X 10^^ G in the Born-Oppenheimer approximation. The variational method is 
employed with a single trial function which includes electronic correlation in the form exp(7ri2), 
where 7 is a variational parameter. It is shown that the quantum numbers of the state of the lowest 
total energy (ground state) depend on the magnetic field strength. The ground state evolves from 
the spin-singlet state for weak magnetic fields B < 5 x 10^ G to a weakly-bound spin-triplet 
state for intermediate fields and, eventually, to a spin- triplet state for 5 x 10^° < ^ < 
4.414 X 10^'^ G. Local stability of the linear parallel configuration with respect to possible small 
deviations is checked. 



*Electronic address: turbiner@nucleares.unam.mx 
^Electronic address: nicolaisOnucleares .unam.mx 
■'•Electronic address: vieyraOnucleares.unam.mx 



1 



I. INTRODUCTION 



The behavior of atoms, molecules and ions placed in a strong magnetic field has attracted 
a significant attention during the last two decades (see, in particular, review papers [1-3]). 
It is motivated by both pure theoretical interest and by possible practical applications in 
astrophysics and solid state physics. In particular, the knowledge of the energy levels can 
be important for interpretation of the spectra of white dwarfs (where a surface magnetic 
field ranges in B 10^ — 10^ G) and neutron stars where a surface magnetic field varies in 
B ^ 10^^ — 10^'^ G, and even can he B ^ 10^'^ — 10^^ G for the case of magnetars. 

Recently, it was announced that in a sufficiently strong magnetic field B > 10^^ G the 
exotic molecular ion Hg"*" can exist in linear configuration with protons situated along the 
magnetic fine [4] (for discussion see a review [3]). In general, it is a metastable long-hving 
system which decays to + p. However, ai B > 10^^ G the ion Hg"^ becomes stable. 
This system does not exist without or for weak magnetic fields. The ion Hg"^ constitutes 
the simplest one-electron polyatomic molecular ion in a strong magnetic field. The Hg^ 
ion has been proposed as being the most abundant chemical compound in the atmosphere 
of the isolated neutron star 1E1207.4-5209 [5]. A detailed review of the current status of 
one-electron molecular systems, both traditional and exotic, that might exist in a magnetic 
field B>10^G can be found in [3]. 

The molecular ion is the simplest stable two-electron polyatomic molecular ion. It 
has a long history since its discovery by J.J. Thomson [6]. Its exceptional importance 
in astrophysics related to interstellar media explains the great interest in this ion from 
astronomy, astrophysics and chemistry communities (for a detailed review, see, [7]). For 
all these reasons, there have been extensive theoretical and experimental works on this 
molecular ion since the pioneering (semi-quantitative) work by Coulson [8]. 

The first variational calculations [9] of the total energy of the molecular ion showed 
that the equilibrium configuration might be either linear or equilateral triangular. However, 
this was not well-established until 1964 [10] when it was shown that the equilibrium con- 
figuration for the state of the lowest total energy is an equilateral triangular configuration, 
while the linear configuration of the Hg ion may occur in excited state(s) . Since that time a 
large number of excited states has been studied [11] (for a general review, see [7]). In partic- 
ular, it has been found that there exists a single spin-triplet state which appears in a linear 
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configuration ^E„. This is also the unique known state of Hg in the hnear configuration. No 
spin-triplet states have been found for a triangular (spacial) configurations so far. 

Although the molecular ion is characterized by the equilateral triangular configuration 
as being the optimal in field-free case, it is expected that in a magnetic field B ^ 0.2 a.u. 
(see below) a linear configuration, parallel to a magnetic field direction, gives the lower total 
energy and becomes the optimal configuration. Somehow, a similar phenomenon already 
happened for the one-electron exotic molecular ion Hg"*" [12] where the optimal configuration 
is triangular at 10^ ^ B < 10^^ G and becomes linear parallel a,t B 10^^ G. It is worth 
noting that for Hg in field-free case the difference between the total energy of the ground state 
(triangular configuration) and of the lowest linear configuration is very small, ~ 0.13 Ry, in 
comparison to characteristic energies in a magnetic field. 

To the best of our knowledge there exists a single attempt to explore the molecular ion 
in a magnetic field [13]. We repeated all numerical calculations of this work following its 
guidelines with use of its formulas (see below. Tables I, V, VI) - in fact, no single number 
from [13] was confirmed. However, in [13] it was made a qualitative statement that with 
a magnetic field increase the transition from equilateral stable equilibrium configuration to 
linear equilibrium configuration may occur. This statement we confirm. Wc predict that 
this transition takes place at a magnetic field ~ 0.2 a.u. A detailed study of a triangular 
configuration and of this transition will be published elsewhere [14]. 

Atomic units are used throughout (h—me—e—l), although energies are expressed in Ryd- 
bergs (Ry). The magnetic field B is given in a.u. with a conversion factor Bq — 2.35 x 10^ G. 

II. GENERALITIES 

Let us consider a system of three protons and two electrons (pppee) placed in a uniform 
constant magnetic field. If for such a system a bound state is developed it corresponds to the 
molecular ion Hg . We assume that the protons are infinitely massive (the Born-Oppenheimer 
approximation of zero order). They are situated along the magnetic field direction forming 
a hnear chain (we call it "the parallel configuration"). The Hamiltonian which describes 
this system when the magnetic field is oriented along the z direction, B = (0, 0, B) is [29] 

k=A,'b,C 
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FIG. 1: Geometrical setting for the ion in linear configuration parallel to a magnetic field 
directed along z-axis. The protons (marked by bullets) are situated on the z-line at distances R± 
from the central proton which is placed at the origin. 

(see Fig. 1 for the geometrical setting and notations), where pi = —iVe is the 3- vector of 
the momentum of the ith electron, the index k runs over protons A, B and C, ri2 is the 
interelectron distance and S = Si + §2 is the operator of the total spin. At is a vector 
potential which corresponds to the constant uniform magnetic field B. It is chosen to be in 
the symmetric gauge, 

^,= ^(Bx r,) = |(-|/,, x,, 0) . (2) 
Finally, the Hamiltonian can be written as 

K=A,B,C 

where — Izi + lz2 ^^^d — + are the z-components of the total angular momentum 
and total spin, respectively, and pt = \/x\^ yj. 

The problem under study is characterized by three conserved quantities: (i) the operator 
of the z-component of the total angular momentum (projection of the angular momentum on 
the magnetic field direction) giving rise to the magnetic quantum number m, (ii) the spatial 
parity operator P{fi — > — ri,f2 — > — r2) which has eigenvalues p — ±l(gerade/ungerade) 
(iii) the operator of the z-component of the total spin (projection of the total spin on the 
magnetic field direction) giving rise to the total spin projection nis- Hence, any eigenstate 



has three exphcit quantum numbers assigned: the magnetic quantum number m, the total 
spin projection nis and the parity p. For the case of two electrons the total spin projection 
nis takes values 0, ±1. 

As a magnetic field increases a contribution from the Zeeman term (interaction of spin 
with magnetic field, B • S) becomes more and more important. It seems natural to assume 
that for small magnetic fields a spin-singlet state is a state of a lowest total energy, while for 
large magnetic fields it should be a spin-triplet state with = — 1, where the electron spins 
are antiparallel to the magnetic field direction B. The total space of eigenstates is spht into 
subspaces (sectors), each of them is characterized by definite values of m, p and m^. It is 
worth noting that the Hamiltonian Ti. is invariant with respect to reflections Zi — > —Zi and 
Z2 — >■ —Z2 (-2-parity operator P^). Hence, any cigenstate is characterized by the quantum 
numbers (Tat = ±1 for positive/negative ^-parity (this symmetry accounts for the interchange 
of the nuclei A and C if they are situated symmetrically with respect to B). 

In order to classify eigenstates we follow the convention widely accepted in molecular 
physics using the quantum numbers m,p and the total spin S without indication to the 
value of rUs. Eventually, the notation is '^^'^^Mp, where 25* -|- 1 is the spin multiplicity which 
is equal to 1 for spin-singlct state {S = 0) and 3 for spin-triplet {S = 1), as for the label 
M we use Greek letters E, 11, A that mark the states with \m\ = 0,1,2,..., respectively, 
and the subscript p (the spatial parity quantum number) takes gerade/ungerade(g'/ii) labels 
describing positive p — +1 and negative p — —1 parity, respectively. There exists a relation 
between the quantum numbers corresponding to the z-parity (interchange of nuclei A and 
C) and the spatial parity: 



Present consideration is limited to the states with magnetic quantum numbers m — 0, —1, —2 
because the total energy of the lowest energy state (the ground state) for any sector with 
m > is always larger than anyone with m < 0. 

As a method to explore the problem we use the variational procedure. The recipe of 
choice of trial functions is based on physical arguments [15]. As a result the trial function 
for a lowest energy state with magnetic quantum number m is chosen in the form 




)m\ im<t>i 'yri2p-oiiriA-a2riB-ci3ric-oi4r2A-ci5r2B-ci6r2C-Bpi^-Bf32^ 



(4) 
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where — ±1 stands for spin singlet (+) and triplet states (— ), while ctjv = 1, —1 stands for 
nuclear gerade and ungerade states, respectively. The P\2 is the permutation operator for 
electrons (1 •(-^ 2) and Pj^, i, j = A, S, C is the operator which interchanges the two protons 
i and j. For S'3-permutationally symmetric case (all protons are identical) Ofq = a^^ = ±1. 
<^i-6! and 7 as well as i?_ are variational parameters. Their total number is eleven. 
It is worth emphasizing that in the trial function (4) the interelectron interaction is included 
explicitly in the exponential form e'^'^^'^ . 

Calculations were performed using the minimization package MINUIT from CERN-LIB. 
Multidimensional integration was carried out using a dynamical partitioning procedure: a 
domain of integration was divided into subdomains following an integrand profile and then 
each subdomain was integrated separately (for details, sec, e.g., [3]). Numerical integration 
was done with a relative accuracy of ~ 10~^ — 10~^ by use of the adaptive DOIFCF routine 
from NAG-LIB. A process of minimization for each given magnetic field and for any partic- 
ular state was quite time-consuming due to a complicated profile of the total energy surface 
in the parameter space but when a minimum is found it took a few minutes to compute a 
variational energy. 
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III. RESULTS 



We carry out a detailed study of E, 11, A low-lying states with a particular emphasis of 
the state which has the lowest total energy for a given magnetic field - the ground state. 

A. m = 

For the case m = we consider four subspaccs in the Hilbert space, 5" = (m^ = 0, spin 
singlet states), S = 1 (spin triplet states) at = —1, = 1 (gerade states) and crjv = — 1 
(ungerade states). 

1. ^Sg state (S = 0, aN = 1) 

For field-free case the system (pppee) in linear configuration (all protons are situated 
on a line, see, Fig.l) the state ^S^ is the lowest total energy state which is characterized 
by a shallow minimum (see, e.g., [7]). However, in spite of developing a minimum for a 
linear chain the system is unstable towards any deviation from linearity. Hence, this state is 
globally unstable. It is worth noting that the true ground state does exist and it corresponds 
to the equilateral triangular configuration (the protons form equilateral triangle) with the 
total energy Ex — — 2.6877Ry and the side of triangle agq — 1.65 a.u. (see [16] and also 
[17])) [30]. A situation is not so different when a magnetic field is not strong, B < 0.2 a.u.: a 
linear parallel configuration with the protons situated along a magnetic line is characterized 
by well-pronounced minimum but a stability towards a deviation from linearity does not 
occur and a global bound state ^E^, does not exist. However, with a magnetic field growth, 
at B > 0.2 a.u. the system (pppee) becomes stable towards small deviations from parallel 
configuration and the ^E^ state exists (being an excited state, see below). 

We made a detailed study of the state ^S^, of the ion in the linear parallel config- 
uration, with a particular emphasis of the symmetric case = R- = R (see Fig. 1), 
as well as small deviations from this configuration in a wide domain of magnetic fields 
< B < 10000 a.u. (see Table I). Finally, for the linear parallel configuration the vari- 
ational trial function ■0*^*"' (4) with ae — 1, ctn — ^ and m — was used. It depends 
on eleven variational parameters. A simple, obvious generalization of (4) is used to study 
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TABLE I: The Hg ion in the state and a comparison with 2e systems H2 and H~: Total Et 
and binding (double-ionization) Ej energies, equihbrium distance Rgq (in a.u.) as well as the total 
energies of final states of dissociation and ionization channels of H3" are shown; all energies are in 
Ry. * the energy and equilibrium distance of Hg" for these magnetic fields is for a case when a linear 
configuration is kept externally (see text), ^ [13], our re-calculations based on the trial function 
from [13] (see text). Total energies for the H2 molecule in state as well as and Hg^ ions 
in lag state in a magnetic field taken from [18], [19] and [20], respectively. Total energies for the 
ground state of the H atom and for the H~ ion in a magnetic field from [21] and [22], respectively. 
The ground state energy of in field- free case from [23] . 
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0.736 
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27.966 


12.034 


0.587 






30.082 








28.317'' 


















28.316^= 




0.576 












100 


177.59 


22.410 


0.336 




181.014 


182.145 


191.361 


190.872 


1000 


1948.41 


51.586 


0.160 






1961.99 


1979.22 


1981.569 


10000 


19891.6 


108.45 


0.083 






19926.25 


19954.60 





"This data can be extracted from [7], p. 427 

slightly deviated configurations when stability of the linear system was checked. 

The variational calculations demonstrate in very clear way the existence of a minimum in 
the total energy surface Et{R+, R-) for the (pppee) system for all magnetic fields ranging 
B — — 10000 a.u. Minimum always corresponds to the symmetric case = = Req of 
the hnear parallel configuration. For B < 0.2 a.u. stabihty is lost with respect to deviations 
from linearity. This indicates a "limited" existence of the molecular ion Hg in the state 
for these magnetic fields. It exists if in some way a linear configuration is supported 
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externally. 

Table I displays the results for the total Et and the double ionization, Ej — 2B — 
Et, energies, as well as for the internuclear equilibrium distance Rgg for the state ^E^,. 
We find that with an increase of the magnetic field strength the total energy grows more 
or less linearly with a magnetic field, the system becomes more and more bound (both 
double ionization and dissociation energies increase) and more compact (the internuclear 
equilibrium distances i?^ and a size of the system L^g — it!+ + i?~ decrease) . 

An important characterization of the system is given by a description of possible dissoci- 
ation and ionization channels together with their behavior as a function of a magnetic field. 
There arc three dominant dissociation channels: (i) Hj]" — H2 + p, (ii) + H and 

(iii) H + p + p (see Table 1) as well as two sub-dominant channels H;^ + p + e 

(ionization) and H3 H + H + p (dissociation) . Last two channels are characterized by 
higher ionization-dissociation energies than the channel H3 — > -l-H and thus they are not 
considered. There are two single-ionization processes — > + p + e and — > Hg"*" + e 
(see Table I). The second one occurs only at S > 10 a. u. where the Rf^ ion can exist, it 
becomes a dominant single-ionization process at B > 10000 a.u. where £^t(H3^) < Et(H^). 
The total energy of the final state compounds after dissociation for different magnetic fields 
is shown in Table I. It is interesting to mention that at B > 100 a.u. the dissociation 
— > Hg"*" -I- e dominates over — > H~ -|- p -|- p. 

A comparison of the total energy of the ground state of H3 for each studied magnetic field 
with the total energy of the products of dissociation or ionization (see Table I) leads to a 
conclusion that the total energy of the ion is always the smallest among them. Thus, the 
B.^ ion in the state is stable for all magnetic fields towards all possible dissociation or 
ionization channels. A smallest dissociation energy corresponds to the channel — > H2-I-P, 
which then is followed by H3 + H. It is worth noting that the largest dissociation 

energy corresponds to the channel — > H~ + p + p. In general, the dissociation energy 
(the difference between the energies of the final and initial states) increases monotonously 
with a magnetic field growth. It is quite interesting that the difference in total energies of 
the final compounds of two major dissociation channels (i) and (ii) grows extremely slow 
with the magnetic field increase reaching 1.1 Ry at B = 100 a.u. 

A conclusion can be drawn that the H3 molecular ion in the state ^12 g exists for B < 
0.2 a.u. if a linear parallel configuration of protons is somehow supported externally, e.g. by 
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placing a system to an (sub)-atomic trap. However, for larger magnetic fields it exists as an 
excited state which is stable towards small deviations from linearity. It is worth noting that 
for the magnetic field B = 0.2 a.u. the total energy well contains at least one longitudinal 

vibrational state. The vibrational energy is calculated following the same procedure which 
is used for Hg"*" ion [20] and it is equal to 0.035 Ry. 

2. state (S = 1, crjv = -IJ 

In field-free case the state of the system (pppee) is (i) the only state of the ion 
in linear configuration which is known so far and also (ii) it is the only known spin-triplet 
state of H3 (for a review of this state see [24] and references therein). For this state several 
vibrational states exist. The linear symmetric configuration i?+ = i?_ is stable towards any 
small deviations, in particular, from linearity. The state is stable with respect to the 
decay H3 -|- H (see [25]). Also there is no decay channel H3 (^E„) H2(^Eg) -|- p. 

A detailed variational study of the ^E^ state of the H3 molecular ion is done for 
< B < 10000 a.u. (see Table II) . It turns out that for all studied magnetic fields the total 
energy surface displays a minimum which corresponds to a linear parallel configuration. Fur- 
thermore, always this minimum appears in the symmetric configuration i?+ = R_ = R. For 
this particular configuration the variational trial function ■z/;*"'^' (4) with de = —1, cttv = — 1 
and m = is used which depends on ten variational parameters. Field-free case is stud- 
ied separately with 23-parametric trial function which is a linear superposition of (4) and 
its three different degenerations [31]. This sufficiently simple function allows to reproduce 
three significant digits in total energy (see Table II). It is separately checked that the linear 
parallel symmetric equilibrium configuration is stable towards all possible small deviations. 

Table II shows the results for the total and the internuclear equilibrium distance Req 
for the ^E^ state for different magnetic fields. With an increase of the magnetic field the 
total energy decreases, the system becomes more bound - double ionization energy increases 
[32] and more compact (the internuclear equilibrium distance decreases). A major emphasis 
of our study of the state ^E^ is the domain 0.2 < B < 20 a.u. where this state becomes 
the ground state of the ion in parallel configuration and likely the global ground state 
of the ion. 

As for the dissociation channel H3 (^E^) H2 (lcrg)-l-H(ls) (with electrons in spin-triplet 
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TABLE II: ion in the state ^T,u' total energy (in Ry), equilibrium distance (in a.u.) and the 
energy of the lowest longitudinal vibrational state E^^'', rotational Eq°* and bending Eq^"-'^. Total 
energy of H^(l(7p) + H(ls) (in Ry) in ground state with spin of each electron antiparallel to B from 
[19] and [21], respectively, shown for comparison. 



B(a.u.) 






rpvib 






rpbend 



Et (at (iTe) + 11(18)) 





-2.2297" 
-2.2322'' 


2.457° 
2.454^ 








-2.2052 


0.1 


-2.3968 


2.416 










0.2 


-2.5991 


2.440 


0.012 


0.0037 


0.014 


-2.5734 


0.5 


-3.0387 


2.273 










1 


-3.6584 


2.125 


0.019 


0.015 


0.028 


-3.6122 


10 


-7.9064 


1.216 


0.048 


0.095 


0.17 


-7.8446 


20 


-10.110 


1.00 


0.063 


0.16 


0.26 


-10.082 


100 


-17.527 


0.645 








-17.855 


1000 


-35.987 


0.372 








-38.01 


10000 


-67.169 


0.235 








-73.75 



°Our calculations (sec text) 
''Rounded data from [25] and [24] 

state) the total energy of the final state is slightly higher than Exiii^) for the magnetic fields 
0.2 ^ B < 20 a.u.; the energy difference varies from 0.03 Ry to 0.06 Ry depending on a 
magnetic field strength, see Table II, remaining very smalL Hence, although Hg (^E^) is 
stable with respect to this dissociation channel it turns out to be a weakly bound state. The 
dissociation may occur at S > 20 a.u. with photon emission at the final state. We do not 
mention a dissociation channel to H2(^S.u) + p due to a probable non-existence of the H2 
molecule in the domain 0.2 < 20 a.u. (see e.g. [18]). 

In the domain 0.2 ^ B < 20 a.u. the total energy well corresponding to the state 
contains at least one longitudinal vibrational state (see Table II). Its energy grows with a 
magnetic field increase. It is calculated the lowest rotational energy as well as the lowest 
bending energies using the same formulas as for Hj"*" [20]. All these energies grow with a 
magnetic field increase. The interesting observation is that for each magnetic field in the 
domain 1 < i? < 20 a.u. the following hierarchy of these energies holds: 

rpvib ^ rnrot ^ rnbend 

contrary to the hierarchy at 0.2 < B < la.u. 
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Hence, the bending energy is the highest to the contrary the hierarchy at the field-free case 
where the longitudinal vibrational energy is the highest (see e.g. [24]), 

rnrot , rnbend - rnvib 
-C'O < -^0 < -^0 ■ 

Finite-proton-mass efi^ects might change the binding energies. So far, it is not completely 
clear how such effects can be calculated quantitatively. At present, a size of their contribution 
might be estimated by values of the energies of the normal modes - the lowest vibrational, 
rotational and bending energies. Their contribution to the binding-dissociation energies 
grows with a magnetic field increase (see Table II) and may reach 10-20% for the magnetic 
fields close to the Schwinger limit (for discussions and references see [3]). 

A comparison of the total energies of the H3 ion for the states ^Eg and (see Tables I 
and II) shows that at S 0.2 a.u. the energy crossing between these two states occurs. 
It implies that for linear parallel configuration the lowest energy state for B < 0.2 a.u. is 
the ^T,g state while for B > 0.2 a.u. the state ^E^ gets the lowest total energy becoming 
the ground state. Hence, one can state that in the domain 0.2 < -B < 20 a.u. the state 
^Ey is the ground state (see below a description of H and A states). In this region of the 
magnetic fields the linear parallel configuration for the state ^E„ is stable towards small 
deviations. It was demonstrated by calculating the corresponding curvatures and then the 
lowest vibrational, rotational and bending energies. It is worth noting that these energies 
(see Table II) turned to be small in comparison to the Et — —Ej which imphes small 
finite-proton-mass effects. 

However, the state as a ground state is weakly bound - energy needed for dissociation 
to H2 (IcTg) -|-H(ls) with electron spins antiparallel to B is very small. This weak boundness 
can be considered as a consequence of the fact that electrons are in the same quantum state, 
thus the Pauli repulsion plays an essential role leading to a large exchange energy. It is 
worth emphasizing that at S ~ 0.2 a.u. the total energy of the global ground state given 
by a triangular configuration [14] coincides approximately to the total energies of the states 
^Eg and ^E„. 
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TABLE III: Hj^ ion in the state ^S^: Total energy (in Ry) and equilibrium distance in (a.u.) (in 
field-free the state does not exist). 



B(a.u.) 






1 


-3.3256 


5.139 


10 


-6.9315 


3.063 


100 


-14.834 


1.958 


1000 


-29.66 


1.35 


10000 


-54.55 


0.94 



3. state (S = 1, un = 1) 

In field-free case the state of the ion in linear configuration does not exist - the 
total energy surface does not reveal a minimum or even irregularity which would correspond 
to this state. However, when a magnetic field is imposed this state may appear. It happens 
already B — 1 a.u. where the total energy surface Et{R+, R-) of this state displays a 
well-pronounced minimum for linear parallel configuration. A detailed variational study of 
the state of the B.^ molecular ion in linear parallel configuration is done for the domain 
1 < B < 10000 a.u. (sec Table III). The trial function V--*™' (4) at ae = -1, cttv = 1 and 
m = is used, it depends on eleven variational parameters. 

The calculations indicate clearly the existence of a minimum in the total energy surface 
Et{R+, R-) of H3 for all studied magnetic fields B — 1 — 10000 a.u. The minimum always 
occurs for the symmetric configuration i?+ — R^ = R. The results are presented in Table III. 
With an increase of the magnetic field strength the total energy decreases. The system 
becomes more bound: the double ionization energy grows [26] . Also the system gradually 
becomes more compact - the internuclear equilibrium distance gradually decreases. 

4. state (S = 0, aN = -l) 

Similar to the state ^T,g in the field-free case the state of the ion in linear 
configuration does not exist. However, when a magnetic field is imposed this state can 
occur. Similar to the state it happens already at -B = 1 a.u. where the total energy 
surface of this state displays a minimum. A detailed variational study of the state ^E„ of 
the H.^ molecular ion in linear parallel configuration is done for 1 < i? < 10000 a.u. (see 
Table IV). The trial function -0*^*"^ (4) at — 1, — —1 and m = is used for it which 
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TABLE IV: ion in the state ^S^: total Et and double- ionization energies Ej (in Ry), and 
equilibrium distance (in a.u.) of H3" (in field-free case this state does not exist). 



B{a.u.) 


Et 


Ei 


Req 


1 


-1.3256 


3.3256 


4.632 


10 


13.0545 


6.9454 


2.563 


100 


185.150 


14.85 


1.651 


1000 


1970.36 


29.64 


1.494 


10000 


19945.6 


54.42 


1.328 



depends on eleven variational parameters. 

The variational calculations indicate clearly the existence of a minimum in the total 
energy surface Et{R+,R-) of H3 for magnetic fields ranging B — 1 — 10000 a.u. The 
minimum always occurs for the symmetric configuration i?+ = R_ = R. In Table IV 
the results for the total Et and double ionization energies {Ej — 2B — Et) as well as 
the internuclear equilibrium distance R^g are shown. With an increase of the magnetic 
field strength the total energy increases, the system becomes more bound (double ionization 
energy increases) and gradually more compact (the internuclear equilibrium distance globally 
decreases) . 

B. m = — 1 

For the case m — —1 four subspaces are studied: S — (spin singlet states) and S — 1 
(spin triplet states) with = —1, and parities — i and — —1, respectively. All 
these states do not exist in the field-free case. 

1. state (S = 1, aN = 1) 

The spin-triplet state ^11^ of the molecular ion in linear configuration does not exist 
for field-free case. However, when a magnetic field is imposed a minimum on the total 
energy surface Et{R+, R-) can occur. This state is studied in the domain of magnetic fields 
10^ < B < 4.414 X 10^^ G using the variational trial function ■0*""' (4) with cXe = — 1, 
aN = 1 and m = —1. It depends on eleven variational parameters. 

The variational calculations indicate clearly the existence of a minimum in the total en- 
ergy surface Et{R+,R-) of for magnetic fields ranging S = 1 - 4.414 x 10^^ G. The 



14 



TABLE V: ion for the state ^n„: total energy Et (in Ry), equilibrium distance Req (in a.u.) and 
the energy of the lowest longitudinal vibrational state E'q**. " [13], Our re-calculations using the 
trial function from [13] (see text). The total energy £'r(H2(^nu)) is from [18] for -8 = 1, 10, 100 a.u., 
while for B = 20, 1000, 10000 a.u. and 4.414 x 10^^ G the total energy is calculated using the present 
technique (it will be described elsewhere). Data for H^(l7ru) and H(ls) from [3] and [26]. 



B(a.u.) 


Ex 


Req 




-BT(H2(3nu)) 


£;T(H+(l7ru)-|-H(ls)) 


ET(H+(lag)-hH(2p_i)) 


1 


-3.036 


1.896 




-2.9686 


-2.6825 


-2.8631 




-2.953 " 














-2.817 


2.040 










5 


-5.654 


1.163 












-5.802 














-5.463 ^ 


1.176 










10 


-7.647 


0.898 




-6.9325 


-6.1980 


-6.5995 




-7.803 " 














-7.307 '' 


0.910 










20 


-9.944 


0.706 


0.135 


-8.934 


-8.036 


-8.582 




-10.475 














-9.752 


0.7 










100 


-18.915 


0.395 


0.343 


-16.473 


-14.452 


-15.547 


1000 


-44.538 


0.183 


1.105 


-35.444 


-31.353 


-33.976 


10000 


-95.214 


0.093 


3.147 


-71.39 


-62.023 


-67.356 


4.414 X 10" G 


-115.19 


0.078 




-84.96 


-73.59 


-79.86 



minimum always corresponds to a linear parallel configuration at i?+ = = R. It was 
investigated its stability towards small deviation from equilibrium in linear parallel config- 
uration. For this state we are not aware how to check quantitatively a stability towards 
deviations from linear parallel configuration. However, physical arguments based on per- 
turbation theory estimates indicate to the stability. Table V contains the results for the 
total Et and the internuclear equilibrium distance Req. With an increase of the magnetic 
field strength the total energy decrease, the system becomes more bound (double ionization 
energy increases [26]) and more compact (the internuclear equilibrium distance decreases). 

The total energy of the final states for the dissociation channels — > H^(^7ru) + H(ls), 
H3 — > H2'(lcrg) -I- H(2p_i) and H3" — > H2(^nu) -|- p with electron spins antiparallel to the 
magnetic field direction for different magnetic fields is shown in Table V. For all studied 
magnetic fields the total energy of both dissociation channels to H^(^7r„) +H(ls), H^(l(Tg) -|- 
H(2p_i) and H2(^n„) are always higher than the total energy of the H3 ion in the 
state. Thus, the ion H^(^n„) is stable towards these decays for all studied magnetic fields. 
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Dominant dissociation channel is Hg — > H2(^nu) +p. For all three channels the dissociation 
energy grows monotonously as a magnetic field increases. For the dominant channel H3 
H2(^nu) + p it reaches 30.3 Ry at the Schwinger hmit 4.414 x 10^^ G, while for the channel 

B.f — > H^(^7ru) + H(ls) for this magnetic field it is required ^ 35 Ry to dissociate. For 
magnetic fields 5x 10^° ^ ^ ^ 4.414x10^^ G there exists at least one longitudinal vibrational 
state (see Table V). 

We made an analysis of the total energies for all spin-triplet states. One can see that 
there is a crossing between the and the states which occurs a,t B 20a.u. It 
shows that the ground state of H3 for B > 20 a.u. is given by the state (see below 
a study of A states which arc characterized by the higher total energies). While the 
state is the ground state for 0.2 < i? < 20 a.u. In Figs. 2 and 3 the evolution of the total 
energy and the equilibrium distance, respectively, of the ground state with the magnetic 
field strength are plotted. The ground state evolves from spin-singlet ^T,g for small magnetic 
fields B < 0.2 a.u. (not shown in Figs. 2,3) to spin-triplet for intermediate fields and 
to spin-triplet state for B > 20 a.u. The total energy decreases monotonously and 
smoothly as magnetic field growth. The equilibrium distance decreases as well, but having 
a discontinuous behavior at B ^ 20 a.u. - in the transition from to states. Similar 
behavior is displayed by < \zi\ >: it reduces monotonously from fa 1.9 a.u. at B — 0.2 a.u. 
to 0.1 a.u. at B — 10000 a.u. with a small discontinuity at B Ki 20 a.u. Perhaps, it is 
worth noting that the average distance between two electrons < ri2 > is also reduced as a 
magnetic field grows in about 20 times between 0.2 a.u. and 10000 a.u. At large magnetic 
fields the transverse size of the electronic cloud coincides approximately with the Larmor 
radius. In Fig. 4 the energy of the lowest longitudinal vibrational state of the ground state 
for 5 X 10^ ^ B < 4.414 x 10^^ G is presented. It grows monotonously as a magnetic 
field increases suffering a discontinuity at S 20 a.u. - in the transition from to 
states. In Figs. 5a,b the valleys and the total energy behavior (profile) along the valley for 
B — 100 a.u. for state, respectively, are shown. Similar behavior takes place for the 
valleys and the total energy profile for ^11^ state for other magnetic fields in the domain 
5 X 10^° < B < AAU X 10^3 Q_ 
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FIG. 2: Ground state energy of viz magnetic field: (stars) and (bullets), a point of 
crossing of these states is marked by empty circle. 
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FIG. 3: Equilibrium distance for the ground state: (stars) and ^11^ (bullets). 
2. state (S = 0, aN = 1) 

A detailed study of the state of the molecular ion in symmetric configuration 
— = R is carried out in the domain of magnetic fields 1 < B < 10000 a.u. (see 
Table VI). The variational trial function ■0*^*"' with — 1, cr^ — 1 and m = — 1 is used for 
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FIG. 4: Energy of the lowest longitudinal vibrational state Eq^^ for the ground state: (stars) 
and (bullets). 

this state, it depends on ten variational parameters. 

The obtained results indicate clearly the existence of a minimum in the total energy Et{R) 
of for all magnetic fields ranging B = 1 — 10000 a.u. Table VI shows the total and 
double ionization energies {Ej — 2B — Et) , as well as the internuclear equilibrium distance 
Req for the state. It is found that with an increase of the magnetic field strength the 
total energy increases, the system becomes more bound (double ionization energies increase) 
and more compact (the internuclear equilibrium distance decreases). 

3. ^Rg state (S = 0, aN = -1) 

It is carried out a detailed study for the state ^Ilg of the molecular ion in symmetric 
configuration — = R in the domain of magnetic fields 1 a.u. < B < 10000 a.u. 
(see Table VII). For this state our variational trial function ■0*""' with Ue — 1, — — 1 
and m — —1 depends on ten variational parameters. The total Et and double ionization 
Ei = 2B — Et energies increase while the internuclear equilibrium distance Rgg decreases as 
a magnetic field grows; the system becomes more bound (double ionization energies increase) 
and more compact (the internuclear equilibrium distance decreases). 
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FIG. 5: Valleys of the total energy (a) and the profile (b) for B = 100 a.u. of the ground state ^n„. 
4. ^Hg state (S = 1, aN = -1) 

A detailed study is carried out for the state '^11^ of the molecular ion in symmetric 
configuration = = Rin the domain of magnetic fields 1 a.u. < B < 10000 a.u. (see 
Table VIII). For this state our variational trial function -0*""' (4) with ae — —1, ctn — — 1 
and m — —1 depends on ten variational parameters. The total Ex energy decreases and 
the double ionization Ej energy [26] increases while the internuclear equilibrium distance 
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TABLE VI 
equilibrium 
(see text). 



ion in the state ^II^j: total Et and double- ionization Ei energies (in Ry) and 
distance Req (in a.u.). " [13], ^ our re-calculations using the trial function from [13] 



B(a.u.) 




Ei 




1 


-0.809 


2.809 


1.995 




-0.561 " 








-0.511 




2.233 


5 


4.747 


5.253 


1.232 




5.025 " 








5.051 ' 




1.262 


10 


13.028 


6.972 


0.967 




13.346 " 








13.346 




0.963 


20 


,30.708 


9.292 


0.750 




31.078 " 








31.081 ' 




0.738 


100 


182.23 


17.77 


0.419 


1000 


1957.77 


42.23 


0.191 


10000 


19909.0 


91.0 


0.098 



TABLE VII: ion in the state ^IIp: total energy Et in Ry and equilibrium distance R^q in a.u. 



B(a.u.) 


Et 


Ei 


Req 


1 


-0.701 


2.701 


3.176 


10 


13.669 


6.331 


1.441 


100 


185.413 


14.587 


0.741 


1000 


1969.42 


30.58 


0.421 


10000 


19942.1 


57.9 


0.273 



Req decreases as a magnetic field grows; the system becomes more bound (double ionization 
energy increases) and more compact (the internuclear equilibrium distance decreases). 



TABLE VIII: ion in the state ^H^: total cikm gy Er in Ry and equilibrium distance Req in a.u. 



B(a.u.) 


Et 




1 


-2.6095 


2.700 


10 


-6.276 


1.487 


100 


-14.429 


0.8.38 


1000 


-30.44 


0.447 


10000 


-57.8 


0.27 
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C. m = -2 



In the m — —2 subspace we study four subspaces: S — (spin singlet states), S — 1 
(spin triplet states with = —1), — 1 and ctjv = —1 and the lowest energy state in 
each of them. All these states do not exist in the field-free case. 

It is carried out a detailed study for the states of the symmetric configuration i?+ = 
R_ = R in the domain of magnetic fields 1 < B < 10000 a.u. For each of these four 
states ^ Ag, ^Ag, ^A^, ^A„ the trial function (4) at m = — 2 depends on ten variational 

parameters. All four states indicate clearly the existence of a minimum in the total energy 
Et{R) of for magnetic fields ranging B — 1 — 10000 a.u. Tables IX - XII show the results. 
For these states with an increase of the magnetic field strength the total energy increases 
for the spin-singlet states and decreases for spin-triplet states, the system becomes more 
bound (double ionization energy increases) and more compact (the internuclear equilibrium 
distance decreases). 

TABLE IX: The ion for the state ^A^: total Et and double- ionization Ej energies in Ry and 
equilibrium distance Reg in a.u. 



B (a.u.) 


Et 


Ei 


Req 


1 


-0.6136 


2.6136 


2.206 


10 


13.499 


6.501 


1.027 


100 


183.325 


16.675 


0.433 


1000 


1960.19 


39.81 


0.191 


10000 


19913.6 


86.4 


0.10 



TABLE X: The Hg" ion in the state ^A^: total energy Et in Ry and equilibrium distance Reg in 
a.u. 



B(a.u.) 




Req 


1 


-2.633 


2.179 


10 


-6.624 


1.013 


100 


-16.92 


0.432 


1000 


-40.38 


0.197 


10000 


-87.49 


0.099 


4.414 X 10" G 


-106.02 


0.09 
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TABLE XI: The Hg ion in the state ^A„: total energy Et in Ry and equihbrium distance Reg in 
a.u. 



B(a.u.) 


Et 


Ei 


Req 


1 


-0.4107 


2.4107 


3.316 


10 


14.281 


5.719 


1.514 


100 


186.602 


13.398 


0.775 


1000 


1972.08 


27.92 


0.401 


10000 


19945.7 


54.3 


0.273 



TABLE XIL The ion in the state ^A^: total energy Et in Ry and equilibrium distance Req in 
a.u. 



B(a.u.) 






1 


-2.443 


4.494 


10 


-5.722 


1.600 


100 


-13.39 


0.804 


1000 


-28.41 


0.449 


10000 


-54.4 


0.28 



IV. CONCLUSION 

We study the low-lying energy states of H3 molecular ion in linear configuration parallel 
to a magnetic field from up to 4.414 x 10^'^ G using the variational method in the Born- 
Oppcnhcimer approximation. The total energy curves display a well pronounced minimum 
at finite internuclear distances at i?+ = R- for the lowest states with magnetic quantum 
numbers m = 0, —1, —2, total spins 5" = 0, 1(777.5 — ~1) and parity p — ±1. A level 
distribution for several magnetic field strengths is shown on Fig. 6. If in field-free case there 
exist only two eigenstates in a linear configuration, but many more states in linear parallel 
configuration can appear when a magnetic field is imposed. 

In general, for all studied states, as the magnetic field increases the equilibrium internu- 
clear distances Req decreases and the system becomes more compact, while the total energies 
of spin-singlet states increase whereas that of spin-triplet states decrease. 

The state of the lowest total energy in linear parallel configuration depends on the mag- 
netic field strength. It evolves from spin-singlet (unstable towards a deviation from linearity) 
for weak magnetic fields B < 0.2 a.u. to spin-triplet (stable towards a deviation from 
linearity) for intermediate fields and eventually to spin-triplet '^n„ state for B > 20 a.u. 
which remains the ground state until the Schwinger limit B — 4.414 x 10^^ G. It is worth 
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emphasizing that for weak magnetic fields, B < 0.2 a.u., the global ground state is given by 
a triangular configuration [14] and then, for larger magnetic fields, the global stable ground 
state corresponds to a hnear parallel configuration [33]. The ion in the state is 
weakly bound. For all studied magnetic fields the total energy surface well corresponding 
to the ground state contains at least one longitudinal vibrational state. 

It is interesting to compare the evolution of the ground state for Hg with magnetic field 
change with that of other two-electron systems (see [27] and references therein). For atomic 
type H~ and He systems there is no domain of magnetic field where the spin-triplet, m — 
state is the ground state: for weak fields the ground state is the spin-singlet, m — state and 
then it becomes the spin-triplet, m = — 1 state for large fields. For the hydrogen molecule 
the state is unbound for all magnetic fields unlike the case of H^. It implies that the 
H2 molecule does not exist as a bound system for 0.18 ^ B < 15.6 a.u., where the unbound 
state has the lowest total energy at infinitely-large distance between protons. A similar 
situation occurs for the He2"^-ion: it does not exist as a bound system for 0.85 ^ B < 
1100 a.u. [28]. 

What is the lowest-lying excited state for weak magnetic fields B < 0.2 a.u. is not clear 
yet. This question, and also the whole domain B < 0.2 a.u., will be studied elsewhere. In 
the domain of magnetic fields 0.2 < B < 5 a.u. the lowest-lying excited state is ^E^, then 
for B > 5 a.u. the lowest- lying excited state is ^n„. For B > 20 a.u., where the ^11^ state 
becomes the ground state, the lowest-lying excited state is ^E„. However, eit B > 1000 a.u. 
until the Schwinger limit the lowest- lying excited state is ^Ag. 

It is interesting to note that at B = 1000 a.u. the Hg ion exists with as the ground 
state {Et = —44.54 a.u.) with two possible excited states: '^A^, {Et = —40.38 a.u.) and 
{Et = —35.99 a.u.) with energies below the threshold of dissociation to H2(^n„) + 
p {Et — —35.44 a.u.). For larger magnetic fields the situation becomes different. For 
instance, at B — 10000 a.u. for the H3 ion {Et — —95.21 a.u.) only one excited state, 
^Ag {Et — —87.45 a.u.), exists with energy below the dissociation threshold to H2(^H„) +p 
{Et = -71.39 a.u.). Similar situation holds up to the Schwinger limit B = 4.414 x 10^^ G: 
a single excited state "^Ag lies below the dissociation threshold. 

It is found that many states in linear configuration which do not exist for B = begin 
to be bound at relatively small magnetic field B 0.2 a.u. A study of the existence of the 
bound states which might appear in a spacial configuration is our goal for a future study. 
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Another goal is related to a study of transition amplitudes for different electronic states. 

Present consideration is based on the use of a simple variational trial function (4) . This 
function can be easily generalized and extended in the same way as was done in a variational 
study of various one-electron systems in a strong magnetic field (see [3]). This will allow 
to improve the present results and might be done in future. However, we are not sure 
that such a study is crucially important. It is related to a fact that typical accuracies in 
astronomical observations of neutron star radiation would not be higher 10~^ — 10~^ unlike 
to spectroscopical accuracies in laboratory where they can be by several orders of magnitude 
higher. 
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FIG. 6: Total energy of the low-lying levels for B = 0,1,10,100 and 1000 a.u. (energy scale is 
kept the same for all presented magnetic fields but the reference points depend on them) 
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